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Dynamical mean-field theory (DMFT) is a non-perturbative technique for the investigation 
of correlated electron systems. Its combination with the local density approximation (LDA) 
has recently led to a material-specific computational scheme for the ab initio investigation 
of correlated electron materials. The set-up of this approach and its application to materials 
such as (Sr,Ca)V03, V2O3, and Cerium is discussed. The calculated spectra are compared 
with the spectroscopically measured electronic excitation spectra. The surprising similarity 
between the spectra of the single-impurity Anderson model and of correlated bulk materials 
is also addressed. 
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1. Introduction 



In many materials with partially filled d and / electron shells, such as the transition 
metals Ti, V, Fe and their oxides or rare-earth metals such as Ce, electrons occupy narrow 
orbitals. This spatial confinement enhances the effect of the Coulomb interaction between the 
electrons, making them "strongly correlated".^ The interplay between the spin, charge and 
orbital degrees of freedom of the correlated d and / electrons and the lattice degrees of free- 
dom leads to a multitude of unusual ordering phenomena at low temperatures. Consequently, 
strongly correlated electron systems are often exceedingly sensitive to small changes in the 
temperature, pressure, magnetic field, doping, and other control parameters. This results, for 
example, in large changes of the resistivity across metal-insulator transitions, of the volume 
across phase transitions, and of the effective electronic masses. Electronic correlations are 
also essential for an understanding of high temperature superconductivity. These properties 
cannot be explained within conventional mean- field, e.g., Hartree-Fock theory, since these 
theories describe the interaction only in an average way and in terms of a static mean field. 
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During the last few years, our understanding of electronic correlation effects has consid- 
erably improved due to the development of dynamical mean-field theory (DMFT);^""'^'' for an 
introduction into DMFT and its applications see Ref."*^^ Within DMFT the electronic lattice 
problem is mapped onto a single- impurity Anderson model with a self-consistency condition. 
This mapping becomes exact in the limit of large coordination number of the lattice^ and al- 
lows one to investigate the dynamics of correlated lattice electrons non-perturbatively at all 
interaction strengths. 

The single-impurity Anderson niodel^^ itself defines one of the fundamental many-body 
problems whose investigation during the last 40 years has led to enormous physical insights 
and important progress in the development of theoretical investigation techniques. ^'^ In par- 
ticular, a new many-body energy scale Tk (the Kondo"*^^ temperature) is known to arise in 
this problem, which is exponentially small in the limit of vanishing hybridization. Below this 
temperature the system can always be understood as a "local Fermi liquid" with strongly 
renormalized quasiparticles.^^ In particular, the single-particle spectrum, i.e., the local den- 
sity of states (DOS), exhibits a generic three peak structure consisting of two broad maxima 
originating from the atomic levels, and a sharp "Kondo" or "Abrikosov-Suhl" resonance at 
the Fermi level, of width Tk- 

We now know that this characteristic three-peak spectrum of the single-impurity model is 
also found in correlated hulk systems, i.e., lattice models and real materials, where the notion 
of a "single impurity" is not applicable. This surprising fact is explained by DMFT. Indeed, 
DMFT is presently the only theoretical approach which can derive the electronic excitation 
spectrum at all energy scales, hence reproducing the incoherent features at high energies 
(Hubbard bands), and the coherent quasiparticle behavior at low energies within the same 
formalism.^^' 

3. Single-Particle Spectrum of Correlated Electrons: One-Band Hubbard Model 

Apparently there is a close relation between quantum impurity physics, with the Kondo 
problem as its paradigm, and the physics of correlated electrons on a lattice as exemplified by 
the Hubbard model. The connection on the technical side is through the mapping of a lattice 
model onto an effective impurity model, as done in DMFT. As a consequence, well-known 
features of the impurity model such as the Kondo resonance reappear in the solution of lattice 
models. In particular, the pinning of the density of states (DOS) at the Fermi level obtained 
within DMFT^ directly corresponds to the Friedel sum-rule for the single-impurity Anderson 
model. This is clearly seen in Fig. 1 where we show the evolution of the local spectral function 
A{ijj) of the DMFT solution for the one-band Hubbard model 
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Fig. 1. Evolution of the T = spectral function of the one-band Hubbard model with a semi-elliptic 
("Bethe") DOS for interaction values U/W = 0,0.2,0.4, . . . , 1.6 {W: band width) calculated with 
the numerical renormalization group. At the critical interaction 11^2 /W ~ 1.47 the metallic so- 
lution disappears and the Mott gap opens. The "Luttinger pinning" at a; = is clearly observable. 



at zero temperature and half filling as a function of local Coulomb repulsion U in units 
of the bandwidth W of non-interacting electrons. Here i denotes the lattice site, t is the 
amplitude for nearest-neighbor hopping on the lattice, and U is the local Hubbard repulsion. 
Magnetic order is assumed to be suppressed ( "frustrated" ) . While at small U the system can 
be described by quasi-particles whose DOS still resembles that of the free electrons, in the 
Mott insulator state the spectrum consists of two separate "Hubbard bands" whose centers 
are separated approximately by the energy U. The latter originate from "atomic" excitations 
at the energies ±[//2 broadened by hopping of electrons away from the atom. At intermediate 
values of U the spectrum then has a characteristic three-peak structure as in the single- 
impurity Anderson model, which includes both the atomic features (i.e., Hubbard bands) and 
the narrow quasi-particle peak at low excitation energies, near cj = O}'^ This corresponds to a 
strongly correlated metal. The structure of the spectrum (lower Hubbard band, quasiparticle 
peak, upper Hubbard band) is quite insensitive to the specific form of the DOS of the non- 
interacting electrons. 

It is important to note that the three-peak spectrum in Fig. 1 originates from a lattice 
model with one type of electrons only. This is in contrast to the single-impurity Anderson 
model whose spectrum shows very similar features, but is due to two types of electrons, namely 
the localized orbital at the impurity site and the free conduction band. Therefore the screening 
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of the magnetic moment which gives rise to the Kondo effect in impurity systems has a different 
origin in lattice systems. Namely, as explained by DMFT, the same sort of electrons provide 
both the local moments and the electrons which screen these moments. DMFT also allows one 
to investigate the periodic Anderson model in which a lattice of localized electrons hybridizes 
with an uncorrclated conduction band (see, for example, Refs.^^). The resulting DOS of the 
localized electrons again shows the typical three-peak structure. There are, however, important 
differences to both the Hubbard model and the single-impurity Anderson model, such as the 
appearance of a hybridization gap at or close to the Fermi level. 

The vanishing of the quasiparticle peak in the Hubbard model signals a "Mott-Hubbard 
metal-insulator transition". This transition between a paramagnetic metal and a paramag- 
netic insulator induced by the Coulomb interaction between the electrons is one of the most 
famous examples of a cooperative phenomenon involving electronic correlations. The question 
concerning the microscopic origin and nature of this transition poses one of the fundamen- 
tal theoretical problems in condensed matter physics.^''' ^-"^ Correlation induced metal- insulator 
transitions (MIT) are found, for example, in transition metal oxides with partially filled bands 
near the Fermi level such as V2O3 doped with Cr^^'^^ (see section 5.2). For these systems band 
theory typically predicts metallic behavior. 

^From a theoretical point of view the Mott transition is a paradigmatic correlation problem 
since it focusses on the competition between kinetic energy and correlation energy of correlated 
electrons in the solid. Here DMFT has led to significant new insights.^'^'^^'^^"^'^ In particular, 
the Mott-transition is found to be of first order at finite temperatures, being associated with 
a hysteresis region in the interaction range Ud < U < Uc2 where Ud and Uc2 ai'c the values 
at which the insulating and metallic solution, respectively, vanishes. The hysteresis region 
terminates at a critical point {U*,T*). For temperatures above T* the transition changes into 
a smooth crossover from a bad metal to a bad insulator. 

The evolution of the spectral function of the half-filled frustrated Hubbard model at finite 
temperatures, T = 0.0276VF, is shown in Fig. 2. This temperature is above the temperature of 
the critical point so that there is no real transition but a crossover from a metallic-like to an 
insulating-like solution. The height of the quasiparticle peak at the Fermi energy is no longer 
fixed at its zero temperature value. This is due to a finite value of the imaginary part of the 
self-energy. The spectral weight of the quasiparticle peak is seen to be gradually redistributed 
and shifted to the upper (lower) edge of the lower (upper) Hubbard band. The inset of Fig. 2 
shows the [/-dependence of the value of the spectral function at zero frequency yl(a; = 0). For 
higher values of U the spectral density at the Fermi level is still finite and vanishes only in 
the limit C/ ^ 00 (or for T ^ 0, provided that U > Uci{T = 0)). 

For the insulating phase DMFT predicts the filling of the Mott-Hubbard gap with in- 
creasing temperature. This is due to the fact that the insulator and the metal are not distinct 
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Fig. 2. Spectral function for the half-filled Hubbard model for various values oiU a.tT = 0.0276W^ > 
T* (in the crossover region). The crossover from the metal to the insulator occurs via a gradual 

suppression of the quasiparticle peak at uj = 0. The inset shows the i7-dependence of A{oj = 0), in 
particular the rapid decrease for U « I.IW [after Bulla, Costi, and VoUhardt^ ''] . 



phases in the crossover regime, implying that the insulator has a finite spectral weight at 
the Fermi level. This behavior has recently been detected experimentally by photoemission 
experiments.^^ 

4. LDA+DMFT 

Although the Hubbard model is able to explain certain basic features of the Mott-Hubbard 
MIT and of the phase diagram of correlated electrons it cannot explain the physics of real 
materials in any detail. Clearly, realistic theories must take into account the explicit electronic 
structure of the systems. 

Recently, the LDA+DMFT method, a new computational scheme that merges electronic 
band structure calculations and the dynamical mean-field theory, was developed. ^^'^^"^^ Start- 
ing from conventional band structure calculations in the local density approximation (LDA) 
the correlations are taken into account by the Hubbard interaction and a Hund's rule coupling 
term. The resulting DMFT equations are solved numerically with a quantum Monte-Carlo 
(QMC) algorithm. By construction, LDA+DMFT includes the correct quasiparticle physics 
and the corresponding energetics and reproduces the LDA results in the limit of weak Coulomb 
interaction U. More importantly, however, LDA+DMFT correctly describes the correlation 
induced dynamics near a Mott-Hubbard MIT and beyond. Thus, LDA+DMFT is able to 
account for the correct physics for all Coulomb interactions and doping levels. 

In the LDA+DMFT approach^^"^^ the LDA band structure is expressed by a one-particle 
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Hamiltonian i^L£)A) ^^'^ is then supplemented by the local Coulomb repulsion U and Hund's 
rule exchange J. This leads to a material specific generalization of the one-band model Hamil- 
tonian (1): 



* i,mj^m' ,(T,a' 

-^LDA = -^LDA - XI XI ^^<^ 

(3) 

i ma 

Here m and m enumerate the three interacting t2g orbitals of the transition metal ion or 
the 4/ orbitals in the case of rare earth elements. The interaction parameters are related by 
V = U — 2J which holds exactly for degenerate orbitals and is a good approximation for the 
t2g- The actual values for U and V can be obtained from an averaged Coulomb parameter 
U and Hund's exchange J, which can be calculated by constrained LDA. The energy term 
containing Ae^ is a shift of the one-particle potential of the interacting orbitals. It cancels the 
Coulomb contribution to the LDA results, and can be calculated by constrained LDA.^° 

Within the LDA-I-DMFT scheme the self-consistency condition connecting the self-energy 
S and the Green function G at frequency uj reads: 

Ggm,g'm'{u;)=^ [ d'k ([ + fxl - H^^j^ik) - E{u;)]-') (4) 

Vb J ^ ^ qm,q'm' 

Here, 1 is the unit matrix, /x the chemical potential, H^^p^{k) is the orbital matrix of 
the LDA Hamiltonian derived, for example, in a linearized muffin-tin orbital (LMTO) basis, 
S(a;) denotes the self-energy matrix which is nonzero only between the interacting orbitals, 
and [...]~^ implies the inversion of the matrix with elements n {=qm), n'{=q'm'), where q 
and m are the indices of the atom in the primitive cell and of the orbital, respectively. The 
integration extends over the Brillouin zone with volume Vb (note that H^^j^ may include 
additional non-interacting orbitals). 

For cubic transition metal oxides Eq. (4) can be simplified to 

if the degenerate t2g orbitals crossing the Fermi level are well separated from the other or- 
bitals.^° For non-cubic V2O3 the degeneracy is lifted. In this case we employ Eq. (5) as an 
approximation, using different Sj„(a;), N^{e) and Gm{ijj) for the three non-degenerate t2g 
orbitals. 

The Hamiltonian (2) is solved within the DMFT using standard quantum Monte-Caxlo 
(QMC) techniques to solve the self-consistency equations. From the imaginary time QMC 
Green function we calculate the physical (real frequency) spectral function with the maximum 
entropy method,^^ using the program by Sandvik. 
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5. Single- Particle Spectrum of Correlated Electrons: Real Materials 

Transition metal oxides are an ideal laboratory for the study of electronic correlations 
in solids. Among these materials, cubic pcrovskites have the simplest crystal structure and 
thus may be viewed as a starting point for understanding the electronic properties of more 
complex systems. Typically, the 3d states in those materials form comparatively narrow bands 
with width W^2 — 3 eV, which leads to strong Coulomb correlations between the electrons. 
Particularly simple are transition metal oxides with a 3d^ configuration since, among others, 
they do not show a complicated multiplet structure. 

Photoemission spectra provide a direct experimental tool to study the electronic structure 
and spectral properties of electronically correlated materials. Intensive experimental investi- 
gations of spectral and transport properties of strongly correlated 3d^ transition metal oxides 
started with investigations by Fujimori et al.'^^ These authors observed a pronounced lower 
Hubbard band in the photoemission spectra (PES) which cannot be explained by conventional 
band structure theory. 

In the following we will employ LDA+DMFT to compute the k-integrated electronic 
spectra of two correlated materials, the 3d^ system (Sr,Ca)V03 and the more complicated 
3d^ system V2O3, using the Hilbert transform of the LDA DOS (see eq. (5)). 

5.1 Sr^Cai-^VOs 

SrVOa and CaVOa are simple transition metal compounds with a 3d^ configuration. The 
main effect of the substitution of Sr ions by the isovalent, but smaller, Ca ions is to decrease 
the V-O-V angle from 9 = 180° in SrVOa to 9 ^ 162° in the orthorhombically distorted 
structure of CaVOa. However, this rather strong bond bending results only in a 4% decrease 
of the one-particle bandwidth W and thus in a correspondingly small increase of the ratio 
U/W as one moves from SrVOa to CaVOs.^^ 

LDA+DMFT(QMC) spectra of SrVOs and CaVOs were calculated by Sekiyama et al.^'^ 
by starting from the respective LDA DOS of the two materials, and are shown in Fig. 3. 
These spectra show genuine correlation effects, i.e., the formation of lower Hubbard bands at 
about 1.5 eV and upper Hubbard bands at about 2.5 eV, with well pronounced quasiparticle 
peaks at the Fermi energy. Therefore both SrVOs and CaVOa are strongly correlated metals. 
The DOS of the two systems shown in Fig. 3 are quite similar. In fact, SrVOs is slightly less 
correlated than CaVOs, in accord with their different LDA bandwidths. The inset of Fig. 3 
shows that the effect of temperature on the spectrum is small for T < 700 K. Spectra of 
SrVOs and CaVOs were also calculated independently by Pavarini et al.^^ The enhancement 
of electronic correlations at the surface of SrVOs as compared to the bulk was studied by 
Liebsch.^^ 

Since the three t2g orbitals of this simple 3d^ material are (almost) degenerate the spectral 
function has the same three-peak structure as that of the one-band Hubbard model shown in 
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Fig. 3. LDA+DMFT(QMC) spectrum of SrVOa (solid line) and CaVOa (dashed line) calculated at 
T=300 K; inset: effect of temperature in the case of CaVOa [after Sekiyama et aZ.^^]. 




Fig. 4. Comparison of the calculated, parameter-free LDA+DMFT(QMC) spectra of SrVOa (solid 
line) and CaVOa (dashed line) with experiment. Left: Bulk-sensitive high-resolution PES (SrVOs: 
circles; CaVOa: rectangles) [after Sekiyama et al.^'^]. Right: Is XAS for SrVOs (diamonds) and 
Cao.9Sro.1VO3 (triangles) .^^ Horizontal line: experimental subtraction of the background intensity. 



Fig. 2. The temperature induced decrease of the quasiparticle peak height is also clearly seen. 
As noted in Sect. 3 the actual form of the spectrum no longer resembles the input (LDA) 
DOS, i.e., it essentially depends only on the first three energy moments of the LDA DOS 
(electron density, average energy, band width). 

In the left panel of Fig. 4 we compare the LDA-FDMFT(QMC) spectra at 300K to ex- 
perimental high-resolution bulk PES. For this purpose we multiplied the theoretical spectra 
with the Fermi function at the experimental temperature (20 K) and Gauss broadened with 
the experimental resolution of 0.1 eV.'^'' The quasiparticle peaks in theory and experiment are 
seen to be in very good agreement. In particular, their height and width are almost identical 
for both SrVOs and CaVOs. The difference in the positions of the lower Hubbard bands may 
be partly due to (i) the subtraction of the (estimated) oxygen contribution which might also 
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remove some 3d spectral weight below —2 eV, and (ii) uncertainties in the ab-initio calculation 
of U. In the right panel of Fig. 4 we compare with XAS data of Inoue et a/.^^'"^^ We consider 
core-hole life time effects by Lorentz broadening the spectrum with 0.2 eV,^^ multiplying with 
the inverse Fermi function (80K), and then Gauss broadening with the experimental resolution 
of 0.36 eV."^^ Again, the overall agreement of the weights and positions of the quasiparticlc 
and upper t2g Hubbard band is good, including the tendencies when going from SrVOs to 
CaVOs (Cao.9Sro.1VO3 in the experiment). For CaVOs the weight of the quasiparticlc peak is 
somewhat lower than in the experiment. In contrast to one-band Hubbard model calculations, 
our material specific results reproduce the strong asymmetry around the Fermi energy w.r.t. 
weights and bandwidths. Our results also give a different interpretation of the XAS than in 
Ref.^^ where the maximum at about 2.5 eV was attributed to an eg band and not to the 
t2g upper Hubbard band. The slight differences in the quasiparticlc peaks (see Fig. 3) lead 
to different effective masses, namely m* /niQ = 2.1 for SrVOs and m* /niQ = 2.4 for CaVOa. 
These theoretical values agree with m*/mo = 2 — 3 for SrVOa and CaVOa as obtained from 
de Haas- van Alphen experiments and thermodynamics.^^' 

5.2 V2O3 

The physical properties of V2O3, and the metal-insulator transition in its paramagnetic 
phase, have been subject of experimental and theoretical studies for more than 30 years. Recent 
advances in experimental PES and the microscopic modeling of correlated electron systems 
by the LDA-I-DMFT approach have led to essential new insights into this correlation-induced 
phenomenon. 

Using the crystal structure of paramagnetic metallic (PM) V2O3 and paramagnetic in- 
sulating (PI) (Vo.962Cro.o38)203, respectively, as input we performed LDA-I-DMFT(QMC) 
calculations with one aig and two degenerate bands. To study the metal-insulator tran- 
sition at experimentally relevant temperatures we calculated at T = TOOK and T = 300 K. 
Since the computational effort is proportional to the low temperature calculations were 
computationally very expensive. Fig. 5 shows the results of our calculations at T = 1160 K, 
T = TOOK, and T = 300K for metallic V2O3, and at T = 1160 K and T = TOOK for insulat- 
ing (Vo.962Cro.o38)203. In the metallic phase, the incoherent features are hardly affected by a 
change in temperature, whereas the quasiparticle peak becomes sharper and thus more pro- 
nounced at lower temperatures. This behavior also occurs in the Anderson impurity model^^ 
and has its origin in the smoothing of the Kondo-Abrikosov-Suhl resonance at temperatures 
larger than the Kondo temperature. However, in contrast to the Anderson impurity model 
this smoothing occurs here at considerably lower temperatures which is apparently an effect 
of the DMFT self-consistency cycle. 

Before focusing on the comparison with experiment, let us briefly discuss a peculiarity of 
the Mott-Hubbard transition which is due to the orbital degrees of freedom.^'^ In the one- 
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(E-Ep) (eV) 

Fig. 5. LDA+DMFT(QMC) spectra for paramagnetic insulating (Vo.962Cro. 038)203 and metallic 
V2O3 at f7 = 5 eV [after Keller et al.^"]. 

band Hubbard model the Mott-Hubbard transition is characterized by the disappearance of 
the quasiparticle weight Z — > 00. However, for inequivalent aig and orbitals as is the case 
considered here only the quasiparticle weight of the orbitals diverges, while that of the aig 
orbital stays finite. Instead, the effective chemical potential at low energies, i.e., /x — in 
Eq. (5), moves outside the non-interacting LDA DOS. As a consequence not the width but 
the height of the aig quasiparticle resonance goes to zero at the Mott-Hubbard transition 
(albeit with a very much reduced width). For the orbitals the transition is characterized 
by a combined shrinking of width and height. In this context we note that the pinning of the 
height of the spectrum at the Fermi energy - valid in DMFT for a single orbital^ (or degenerate 
orbitals) - does not hold in the case of inequivalent orbitals. In our case, the volume enclosed 
by the aig orbital shrinks at the expense of the orbitals. Only the total volume remains 
constant in accord with Luttinger's theorem. 

To be able to compare with experimental photoemission spectra, the LDA-I-DMFT re- 
sults^^ were multiplied with the Fermi function at the experimental temperature (T ~ 180 K) 
and broadened with a 0.09 eV Gaussian to account for the experimental resolution. The same 
procedure was used for the comparison with x-ray spectroscopy data (with an inverse Fermi 
function at T = 300 K and a broadening of 0.2 eV taken from experiment). On the experimen- 
tal side, the PES of Refs.^^'^'^ were corrected for the inelastic Shirley-type background which 
also removes the 0-2p contribution. All experimental and theoretical curves were normalized 
to yield the same area (which is a measure of the occupation of the vanadium t2g bands). 

In Fig. 6, the LDA-I-DMFT results at 300 K are compared with early photoemission spectra 
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Fig. 6. Comparison of LDA+DMFT(QMC) results at T = 300 K with photoemission data by 
Schramme et al.^'^ and Mo et al.'^^ for metallic V2O3 [after Keller et al.'^'^]. 

by Schramme^'^ and recent high-resolution bulk-sensitive photoemission spectra by Mo et al.'^^ 
The strong difference between the experimental results is now known to be due to the distinct 
surface sensitivity of the earlier data. In fact, the photoemission data by Mo et al."^^ obtained 
at hv = 700 eV and T = 175 K exhibit, for the first time, a pronounced quasiparticle peak. 
This is in good qualitative agreement with our low temperature calculations. However, the 
experimental quasiparticle peak has more spectral weight. The origin for this discrepancy, for 
a system as close to a Mott transition as V2O3, is presently not clear. 

While the comparison with PES data provides important insight into the physics of V2O3, 
more than half of the theoretical spectrum lies above Ep- For this region we compare our 
results at 1160 K, 700 K, and 300 K with O Is X-ray absorption spectra (XAS)^^ for V2O3 at 
300 K by Miiller et al.^^ (see Fig. 7). The theoretical spectra above Ep are found to be almost 
independent of temperature. Namely, there is a shoulder at higher temperatures developing 
into a small peak at low temperatures (300 K)) which is the residue of the quasiparticle peak. 
Furthermore, at 1.1 eV there is a rather narrow peak, and at about 4.2 eV a broad peak. These 
two peaks are parts of the upper Hubbard band, and are due to the presence of more than 
one type of correlated electron in the problem {aig, eg) with Hund's rule coupling J. The 
latter leads to a splitting of the upper Hubbard band. Hence, the relative position of those 
two peaks depends sensitively on the value of J. In particular, a slightly smaller value of J 
will make the agreement with experiment even better. 

6. Volume collapse in Cerium 

Cerium (Ce) exhibits a transition from the 7- to the a-phase with increasing pressure or 
decreasing temperature. This transition is accompanied by an unusually large volume change 
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of 15%,^^ much larger than the 1-2% volume change in V2O3. The 7-phase may also be 
prepared in metastable form at room temperature in which case the 7-a transition occurs 
under pressure at this temperature.^^ Similar volume collapse transitions are observed under 
pressure in Pr and Gd (for a recent review, see McMahan et a/.^^). It is generally believed that 
these transitions arise from changes in the degree of 4/ electron correlations, as is reflected in 
both the Mott transition^^ and the Kondo volume collapse (KVC)^^'^'' models. 

The Mott transition model envisions a change from itinerant, bonding character of the 4/- 
electrons in the a-phase to non-bonding, localized character in the 7-phase, driven by changes 
in the 4/-4/ inter-site hybridization. Thus, as the ratio of the 4/ Coulomb interaction to the 
4/-bandwidth increases, a Mott transition occurs to the 7-phase, similar to the Mott-Hubbard 
transition of the 3d-electrons in V2O3. 

The Kondo volume collapse^^ scenario ascribes the collapse to a strong change in the 
energy scale associated with the screening of the local 4/-moment by conduction electrons 
(Kondo screening), which is accompanied by the appearance of an Kondo- Abrikosov-Suhl- 
like quasiparticle peak at the Fermi level. Indeed, the 7-a-transition can be described using 
a single- impurity model where the model parameters are determined from DFT/LDA and 
spectroscopy.^"^ 

In the KVC model the change of the Ce-4/-electron spectrum across the transition is, in 
principle, very similar to that in the Mott scenario, i.e., there will be a strong reduction of 
the spectral weight at the Fermi when going from the a- to the 7-phase. The subtle difference 
comes about by the 7-phase having metallic /-spectra with a strongly enhanced effective mass 
as in a heavy fermion system, in contrast to the /-spectra characteristic of an insulator in 
the case of the Mott scenario. The /-spectra in the Kondo picture also exhibit Hubbard side- 
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bands not only in the 7-phase, but in the ct-phase as well, at least close to the transition. 
While local-density and static mean-field theories used in the Mott transition model up to 
now^^ correctly yield the Fermi-level peaks in the /-spectra for the a-phase, they cannot 
reproduce the Hubbard side-bands since these treatments neglect genuine correlation effects. 
By contrast, DMFT solutions of both Hubbard and periodic Anderson models do exhibit the 
Hubbard side-bands in the a-like regimes.®^ 

Typically, the Hubbard model and the periodic Anderson model are considered paradigms 
for the Mott and KVC model, respectively. Although both models describe quite different 
physical situations it was recently shown that they lead to a surprisingly similar behavior 
at finite temperatures. Namely, for increasing Coulomb interaction the spectrum and local 
magnetic moment show very similar features. This is also the case for the phase diagram and 
the charge compressibility of the periodic Anderson model with nearest neighbor hybridiza- 
tion.^^' ^'^ From this point of view the two scenarios are no longer really distinct, at least at 
temperatures relevant for the description of the a-7 transition. 

For a realistic calculation of the Ce a-7 transition, we employ the full Hamiltonian calcu- 
lation described in Sect. 4 where the parameters entering the one-particle Hamiltonian were 
calculated by LDA and the 4/ Coulomb interaction U along with the associated 4/ site energy 
shift by constrained LDA (for details of the two calculations presented here see Refs.^^'^'^'^-'" 
and Ref.59). 

The LDA-I-DMFT(QMC) spectral evolution of the Cc 4/-elcctrons is presented in the left 
panel of Fig. 8. At a volume per atom V = 20 A^, this figure shows that almost the entire 
spectral weight lies in a large quasiparticlc peak with a center of gravity slightly above the 
chemical potential. This is similar to the LDA solution; however, a weak upper Hubbard band 
is also present even at this small volume. At the volumes 29 and 34 A'^ which approximately 
bracket the a-7 transition, the spectrum has a three peak structure. Finally, at V = 46 A"^ the 
central peak has disappeared, leaving only the lower and upper Hubbard bands. An important 
difference to V2O3 is the metallic feature of the spd-spectrum. Thus Ce remains a metal across 
the transition which is monitored by a vanishing 4/ quasiparticle resonance. 

To study the energetic changes associated with the rapid change of the quasiparticle weight 
at the Fermi energy, we calculate the DMFT energy per site for the model Hamiltonian (3) 

Ebmft = ^ j;TV(ifODA(k)Gk(i^„))e^'^»°^ + C//d. (6) 

nkiT 

Here, Tr denotes the trace over the 16 x 16 matrices, T the temperature, N the number of k 
points, Gk the Green function matrix w.r.t. the orbital indices, i?Lj-,^(k) the LDA one-particle 
matrix , and 

is a generalization of the one-band double occupation for multi-band models. The prime on 
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Fig. 8. Spectral functions for Ce as calculated by LDA-I-DMFT(QMC). Left: 4/ spectral function 

A{uj) at different volumes and T = 632K {uj = 0: chemical potential; curves are offset as indicated; 
Ar = 0.11cV^^). Right: Total LDA+DMFT spcZf-spcctrum (solid line) shown in comparison with 
the measured PES^'' and BIS^^ (circles) for a- (upper part) and 7-Ce (lower part) at T = 580 K 
[after McMahan et al.'^^]. 



the sum indicates that at least two of the indices of an operator have to be different. 

Fig. 9a shows our calculated DMFT(QMC) energies -Edmft as a function of atomic vol- 
ume at three temperatures relative to the paramagnetic Hartree-Fock (HF) energies -Epmhf 
of the Hamiltonian (3), i.e., the energy contribution due to electronic correlations. We also 
present the HF energies of a polarized solution which basically represent a non-self-consistent 
LDA-)-U calculation and reproduce -Edmft at large volumes and low temperatures. With de- 
creasing volume the DMFT energies bend away from the polarized HF solutions. Thus, at 
r = 0.054 eV;^ 600 K a region of negative curvature in -Edmet — -E^pmhf is evident within the 
observed two phase region (arrows). 

Fig. 9b presents the calculated LDA-I-DMFT total energy Etot{T) = Ei,BA{T)-\-EBMFT{T)- 
-E'mLDA(T') where -EmLDA is the energy of an LDA-like solution of the Hamiltonian (3).^^ 
Since both E'lda and -Epmhe — ^mLDA have positive curvature throughout the volume range 
considered, it is the negative curvature of the correlation energy in Fig. 9a which leads to 
the dramatic depression of the LDA-I-DMFT total energies in the range F =26-28 for 
decreasing temperature, in contrast to the smaller changes near V=34A^ in Fig. 9b. This 
trend is consistent with a double-well structure emerging at still lower temperatures (which is 
prohibitively expensive for QMC simulations) , and with it a first-order volume collapse. This 
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Fig. 9. (a) Correlation energy -Bdmft— -Epmhf shown as a function of atomic volume (symbols) and 
Hartree-Fock (HF) energy -Eafhf— -Epmhf of a polarized solution (dotted lines); arrows: observed 
volume collapse from the a- to the 7-phase. The correlation energy is seen to bend away from the 
HF energy of the polarized solution in the region of the transition, (b) Total energy near V = 26-28 
A3 ( same symbols as in (a)); curves at T = 0.544 eV were shifted downwards by — 0.5eV to match 
the energy range [after Held et a/.^°]. 



is in reasonable agreement with the experimentally found volume collapse. Other physical 
quantities like the free energy and the specific heat are discussed in Ref.^^ We note that a 
similar scenario has recently been proposed by Savrasov et al.^^ for the S-a transition in Pu 
on the basis of LDA+DMFT calculations. These authors solve the DMFT equations by an 
ansatz inspired by iterated perturbation theory ^'^ and include a modification of the DFT/LDA 
step to account for the density changes introduced by the DMFT. 

The comparison of various static physical quantities of Ce calculated by 
LDA+DMFT(QMC) and LDA+DMFT(NCA)59 ^j^j^ experiment shows fair to good 
agreement in the overall behavior and, except for the susceptibility, even in absolute values. 

Finally, in the right panel of Fig. 8 we compare the spectral functions calculated by 
LDA+DMFT (QMC) (see left panel of Fig. 8) with experiment.^^ The main contribution to 
the PES of a-Ce (right panel of Fig. 8) is seen to come from the energy range between 
3 eV and 7 eV, which is attributed to 4/^ final state multiplets. In the calculated spectrum 
all excitations to 4/^ states are described by the featureless upper Hubbard band. As a 
consequence of the simplified interaction model all doubly occupied states are degenerate. 
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This shortcoming in our calculation is responsible for the sharply peaked main structure. The 
neglected exchange interaction would produce a multiplet structure closer to experiment. The 
calculated /-spectrum shows a sharp quasiparticle resonance slightly above the Fermi energy, 
which is the result of the formation of a singlet state between /- and conduction states. We 
thus suggest that the spectral weight seen in the experiment is a result of this quasiparticle 
resonance. In the lower part of the right panel of Fig. 8, a comparison between experiment and 
our calculation for 7-Cc is shown. The most striking difference to the result for a-Ce (upper 
part of the right panel of Fig. 8) is the absence of the quasiparticle resonance in the 7-phase 
which is in agreement with our calculations. Nonetheless 7-Cc remains metallic with spectral 
weight arising from the spd-electrons at the Fermi energy. Altogether, one can say that the 
agreement with the experimental spectrum is very good — comparable to the accuracy of 
LDA for much simpler systems. 

7. Conclusion 

In this paper we discussed the set-up and several applications of the computational scheme 
LDA+DMFT which merges two non-perturbative, complementary techniques for the theo- 
retical investigation of many-body systems in solid state physics. Using the band structure 
calculated within local density approximation (LDA) as input, the missing electronic correla- 
tions are introduced by dynamical mean-field theory (DMFT). Thereby LDA-I-DMFT allows 
one to perform ab initio calculations of real materials with strongly correlated electrons, 
i.e., electronic systems close to a Mott-Hubbard metal-insulator transition, heavy fermions, 
and /-electron materials. The physical properties of such systems are characterized by the 
correlation-induced generation of small, Kondo-like energy scales which are missing in the 
LDA and which make the application of genuine many-body techniques necessary. 

On a technical level LDA-I-DMFT requires the solution of an effective self-consistent, 
multi-band Anderson impurity problem by numerical methods, e.g., QMC. The investigation 
of quantum impurity problems in the last 40 years, in particular the development of theoretical 
and numerical methods to solve the Kondo problem, is therefore a prerequisite for a successful 
treatment of lattice problems within DMFT. 

The application of LDA-I-DMFT to the transition metal oxides Sr-rCai_a;V03, V2O3 as 
well as elemental Cerium discussed in this paper yield, for example, spectral functions which 
can be compared with spectroscopic measurements. Remarkably good agreement was found 
with photoemission and x-ray absorption data of bulk-sensitive experiments. The spectral 
function of correlated metals are quite generally characterized by a lower Hubbard band, a 
quasiparticle peak near the Fermi energy, and an upper Hubbard band; the latter is split by 
Hund's rule coupling in the case of more than one type of orbital. In the upper Hubbard 
band the orbital structure of a correlated material is therefore particulary clearly displayed. 
This spectrum is similar to that of a localized orbital at the impurity site hybridizing with 
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free conduction electrons ("Kondo problem") although its physical origin is quite different. 
Namely, in correlated bulk materials the same sort of electrons can provide both the local 
moments and their screening. DMFT is able to give a consistent explanation of this surprising 
feature. 

In spite of the remarkable successes of the LDA+DMFT approach for the investigation 
of correlation effects in real materials the method still needs to be improved. Namely, at 
present the LDA band structure serves only as input information for the DMFT, but there 
is no feedback from DMFT to LDA. Since correlation effects can, in principle, change the 
charge distribution on which the LDA band structure depends one needs to feed the changes 
computed by DMFT back into LDA, and repeat the calculation until convergence is reached 
in both parts; for a first implementation see Ref.^® 

Realistic calculations should not only include the orbitals of the correlated electrons, 
but all orbitals. While this is already the case in the Cerium calculations presented here, it is 
difficult to do in the case of transition metal oxides. To overcome these difficulties, an extended 
computational scheme has recently been developed in Wannier basis. It will be applied in 
our future investigations to perform full orbital calculations for correlated materials. 

Besides improving the self-consistency of the LDA+DMFT method, there are also at- 
tempts to improve on the LDA and the DMFT part. For example, instead of LDA the so- 
called GW approximation can be employed. The main advantage is that GW is a purely 
diagrammatic approach. Combined with DMFT, GW+DMFT includes the full contribution 
of the Hartrec diagram, the Fock diagram, and the bubble diagrams for the screening of the 
Coulomb interaction as well as the (DMFT) local contribution of all Feynman diagrams. A 
first simplified implementation was reported by Bicrmann et al.^^ 

Furthermore, non-local correlations may be taken into account by cluster extensions of 
DMFT69,70 (fQj. 

a review sec Maier et al7^), where instead of one Anderson impurity one now 
has several coupled sites hybridizing with the fermionic bath. 
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